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In adiabatic Cooper pair pumps, operated by means of gate voltage modulation only, the quantiza- 
tion of the pumped charge during a cycle is limited due to the quantum coherence of the macroscopic 
superconducting wave function. In this work we show that it is possible to obtain very accurate 
pumps in the non-adiabatic regime by a suitable choice of the shape of the gate voltage pulses. We 
determine the shape of these pulses by applying quantum optimal control theory to this problem. 
In the optimal case the error, with respect to the quantized value, can be as small as of the or- 
der of 10"^e: the error is reduced by up to five orders of magnitude with respect to the adiabatic 
pumping. In order to test the experimental feasibility of this approach we consider the effect of 
charge noise and the deformations of the optimal pulse shapes on the accuracy of the pump. Charge 
noise is assumed to be induced by random background charges in the substrate, responsible for the 
observed 1// noise. Inaccuracies in the pulse shaping are described by assuming a finite bandwidth 
for the pulse generator. In realistic cases the error increases at most of one order of magnitude as 
compared to the optimal case. Our results are promising for the realization of accurate and fast 
superconducting pumps. 

PACS numbers: 73.23.-b, 74.50.+r, 74.78.Na 
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I. INTRODUCTION 



A dc current can be generated in a mesoscopic circuit 
connected to two electrodes even in the absence of a bias 
voltage if a set of external parameters, e.g. gate volt- 
ages, are changed periodically in time^. This mechanism 
is known as charge pumping. If the external parame- 
ters are changed slowly as compared to the typical time 
scales of the system the pumping is adiabatic. In this 
case the charge pumped after one cycle does not depend 
on the detailed timing of the cycle, but only on its ge- 
ometrical properties. Charge pumping can be realized 
in a variety of different situations. In systems consist- 
ing of mesoscopic phase-coherent conductors parametric 
pumping can be achieved through a periodic modulation 
of the phases of the scattering matrix. This regime has 
been studied extensively both for metallic conductors, see 
e.g. Refs. [3Hi4,v.5,i6,i7. and hybrid systems*"^"^° contain- 
ing superconducting terminals. In the opposite regime of 
systems consisting of quantum dots connected through 
tunnel junctions, charge pumping is achieved by the peri- 
odic modulation of the Coulomb blockadeii. In this case 
phase coherence is irrelevant and the number of electrons 
transfered per cycle is approximately quantized: the gen- 
erated current / is related to the frequency of cycle / via 
/ ~ (ne)f, where n is the number of electrons transfered 
in each cycle and e is the electron charge. Experimental 
evidence for parametric charge pumping in normal metal- 
lic systems has been demonstrated for the first time in 
Refs. [mi 

The situation is radically different if superconducting 
islands are considered. Here at low temperatures pump- 



ing is due to the transport of Cooper pairs^^ii^. A Cooper 
pair pump can be realized by an array of Josephson junc- 
tions connected to two superconducting reservoirs, kept 
at a fixed phase bias (p^. In this case, even in those situ- 
ations where pumping is associated to a periodic modula- 
tion of the Coulomb blockade, superconducting phase co- 
herence is fundamentaliiii^. Moreover, in addition to the 
dependence of the pumped charge on the characteristics 
of the cycle, in superconducting pumps there is a depen- 
dence on the superconducting phase difference between 
the reservoirs. The geometric nature of the pumped 
charge has been analyzed both in the Abeliani^i2£i2ii^ 
and non-Abelian23, cases thus opening the possibility to 
experimentally detect geometric phases in superconduct- 
ing nanocircuitg- ^'^^ . Indeed an experimental detection 
of the Berry phase by means of Cooper pair pumping has 
recently been reported^^. 

In addition to the importance of addressing fundamen- 
tal questions related to the quantum mechanical behav- 
ior of macroscopic systems, accurate charge pumping can 
be used for metrological purposes. In the case of single- 
electron pumps the transfered charge, at frequency / of 
a few MHz, has reached such an accuracy (uncertainty 
of 10~*) to make a new metrological standard of capac- 
itance possible^ i^^'^^ . By pumping a certain number of 
electrons onto a capacitor, and measuring the resulting 
voltage, one can measure capacities of the order of IpF. 
On the contrary, the frequencies, and consequently the 
current intensities, in single-electron pumps has not been 
sufficiently high for setting a standard of current. A su- 
perconducting Cooper pair pump, in principle, would al- 
low for higher frequencies although several effects such as 
Landau-Zener tunneling, supercurrent leakage through 
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FIG. 1: Schematic drawing of a Cooper pair pump consisting 
of three Josephson junctions and two gated Cooper pair boxes. 
Ci and Eji are the capacitances and Josephson energies of the 
junctions. Cgi and Vgi are gate capacitances and voltages, 
is the overall phase difference. By a proper modulation of the 
gate voltages a given charge is moved from, say, the left to the 
right electrode. In general the pumped charge will depend on 
the modulation shape and on the phase difference ip. In the 
protocol considered here the dependence on tp is absent. 



the pump, and coherent corrections (which are of cru- 
cial importance to reveal geometric phases) would lead to 
all sorts of inaccuracies in Cooper pair pumping. Some 
works have been done to achieve more accurate pump- 
ing in adiabatic regime by optimizing the design of the 
pumpi^i22ii2£. In this paper we follow a different ap- 
proach. We will design very accurate Cooper pair pumps 
by optimizing the pulse shapes of the gate voltages during 
the cycle. Most importantly we are not bound to work 
in the adiabatic limit therefore increasing, at the same 
time, both the accuracy and the frequency of operation. 
The theoretical framework which will apply to achieve 
this goal is that of optimal control theor y^^'^^i'^^'^'^ . This 
approach was successfully applied to superconducting 
qubits to improve one- and two-qubit gates-'':-^ and it 



was shown that it can help in reducing the error in the 
gate operation by several order of magnitudes. As we 
will describe in the rest of the paper, optimal control is 
also helpful in the case of Cooper pair pumping where, 
for the case we consider, the error is reduced up to five 
orders of magnitude. 

The paper is organized as follows. In the next Section 
we will describe the model for the Cooper pair pump 
used in the rest of the paper. We will then give (Sec. lIIIj) 
a brief introduction to the optimal control algorithm 
employed in this work. In Sec. IIVI we present the results 
obtained for our non adiabatic optimal pump. After 
having discussed the achieved accuracy (Sec. IIV Ap we 
analyze in details the effect of possible imperfections 
in the pulse shapes (Sec. lIVBp . the effect of noise 
(Sec. IIV Cp . as well as the effect of unavoidable varia- 
tions in charging and Josephson energy from junction to 
junction (Sec. lIVPp . in order to test the robustness of 
the pump. The last Section is devoted to the conclusions 
and to a comparison with other proposals to realize 
accurate pumps. 

II. THE COOPER PAIR PUMP 

The three-junction Cooper pair pump considered in 
the paper is shown in Fig. [1] It consists of an array of 
three Josephson junctions with two gated islands; Eji 
and Ci are, respectively, Josephson energy and capacity 
of ith junction, while Cgi and Vgi are, respectively, gate 
capacitance and gate voltage applied on the ith island. 
The number of extra Cooper pairs on the zth island is 
denoted with n,; and ip is the overall phase difference be- 
tween the two superconducting electrodes. In the case 
of a uniform array, which we consider here, the Hamil- 
tonian, in the basis of charge eigenstates | n), is given 



H 



[(ni-qi)'^ + {n2-q2f + ("i - gi) (ns - gs)] | n){n \ - _i^(e^f | n + 4><n \+h.c.) \ (1) 
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k=l 



In this equation n= (711,^2) and q= (91,(72) specify, 
respectively, the number of Cooper pairs on each island 
and the normalized gate charges [qu = —Cg^kVg^k/C^e)). 
Tunneling of one Cooper pair through the fcth junction 
changes the number of pairs on fcth and (fc — l)th island, 

so that the only non-zero components of 5k are {Sk)k = 1 

and {Sk)k-i = —1- Moreover the forward (backward) 
tunneling of one pair through each junction is related 
to a ^ {—j) phase difference. The pump operates in the 
charging limit, i.e. Ej ^ Ec- As discussed in Ref. II Tithe 
Cooper pair pump is operated by changing periodically 



in time the two gate charges qi and q2. 

It is useful to first analyze the case in which the pump 
is driven adiabatically. In this case the time dependence 
of the gate charge is not important, what matters is the 
path which is followed in the parameter space. As was 
discussed extensively in the past (see e.g. Ref. [l^flTt ) 
it is important that the cycle encloses the triple point 
degeneracy (qi = 52 = 1/3) in the stability diagram as 
described in Fig[2l The various regions in Fig[2] indicate 
the corresponding ground state of the charging part of 
the Hamiltonian. After a cycle, for example the one in- 
dicated in the figure, the state of the pump goes back 
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FIG. 2: (Color on line) The stability diagram of a three- 
junction Cooper pair pump in the {qi,q2) space. Each 
hexagon corresponds to a given charge (eigen)state (for Ej — 
0) of the system (711,712). In the adiabatic regime charge is 
pumped when the gate charges are varied along, for example, 
the circular path with radius qo centered at the degeneracy 
point (gi = q2 = 1/3). 



to its initial situation but (approximately) one Cooper 
pair has been transferred through the system. As it has 
been discussed ir^ii, this pump does not lead to a suffi- 
cient quantization accuracy: the error scales as Ej/Ec- 
A smaller error can be achieved either by increasing the 
number of junctions or by pumping by means of gate volt- 



age and flux modulation as in the Cooper pair sluicei^. 

In the present work we take a different approach to 
optimize the pump. Our proposed device operates in 
the non-adiabatic limit, in which the time dependence 
of the pulses is important. In order to have an accurate 
pump wc will then optimize the pulse shapes by means 
of quantum optimal control (see next Section). These 
ingredients allow to construct a fast and accurate Cooper 
pair pump. 

As we will show in the next section, in order to use the 
quantum optimal control, one needs to have a desired fi- 
nal state differing from the initial state and at least one 
parameter as a control in the Hamiltonian of the sys- 
tem. Since in pumping the initial and final states of the 
system are the same apart from a phase we introduce a 
new quantum number, the counter \m). This counter is 
a passive detector which clicks each time a Cooper pair 
goes through the junction where the counter is set. In 
the following we assume that the counter is placed in 
the last junction of the array, therefore if one pair passes 
this junction from left (right) to right (left) the counter 
will change by 1 (—1). We emphasize that the counter 
considered is merely a fictitious element which we use to 
compute the pumped charge. No back action can be ex- 
pected since we are not describing a physical detector. If 
one denotes the state of the array by |<I'(t)) and the state 
of the counter by |m), where m = 0, ±1, ±2..., the state 
of the system will be |$(f)) ® |?7i). In the Hilbert space of 
charge states plus the counter index the Hamiltonian of 
a three-junction uniform pump has the following form: 



J 



TO 



^Ej ^ |^e*3 {\ni + 1,712, to) + |"i - 1,"2 + 1,to) -I- \ni,n2 -l,m+ 1)) (rii, 



n2,m\ 



(2) 



As it is clear in the Hamiltonian ([2]), one can change 
the state of the counter without changing the energy so 
that all states |<i>(t)) (g) |to) with different values of m are 
degenerate. 

III. OPTIMAL QUANTUM CONTROL 

Quantum optimal control algorithms^i^^i^i^i are de- 
signed to lead a quantum system with state from 
an initial state IV'(O)) = \tpini) to a target final state 
hA/in) at time T with a good fidelity where J- = 
\{'>Pfin\'4'{T))\'^. The algorithm employed in this work 
is called immediate feedback control and is guaranteed 
to give a fidelity improvement at each iteration^l. The 
procedure is as follows: Assume that {uj{t)} is the set 



of control parameters which the Hamiltonian depends on. 
First the state of the system \ip{t)) is evolved in time with 
the initial condition \ip{0)) — \^ini) and an initial guess 
{u^^\t)} for control parameters, giving rise to \'>p{T)) af- 
ter time T. At this point an iterative algorithm starts, 
aiming at improving the fidelity by adding a correction 
to control parameters in each step. In the nth step of 
this iterative algorithm 

• The auxiliary state \xiT)) = \->pfin){ipfin\'ip{T)) is 
evolved backwards in time reaching |x(0)). IxC^) 
can be interpreted as the part of the desired final 
state \ipfin) which has been reached. 

• The states |x(0)) and \ip{0)) are evolved forward 
in time with control parameters {u^"^\t)} and 
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{u^"'~^^\t)} , respectively. Here, 



,(«+!) 



,(«) 



{xit)\ 



dH 

duj {t) 



(3) 



are updated control parameters and X(t) is a weight 
function used to fix initial and final conditions on 
the control parameters. 

These two steps are repeated until the desired value of 
fidelity is obtained. 

In this work we will use the described quantum op- 
timal control algorithm to drive the system from the 
initial state \'^ini) = |G(0)) |0) to the final state 
|*/i„) = |G(0)) (g) using the Hamihonian ^ for 
the time evolution where \G{0)) is the ground state of 
the Hamiltonian ([T|) at time t = 0. The initial state 
\G{0)) (g) |m = 0) can be achieved by suddenly couphng 
the array to the electrodes. The normalized gate charges 
qi and (72 will be the control parameters. In all of the 
numerical simulations we performed we assume ip = so 
that supercurrent is zero. Moreover, in the adiabatic case 
and if = the errors are most severe, so that our settings 
describe the worst case scenario to test optimal control 
theory. The expectation value of the counter operator 
m = X^m"^ \m){m\ represents the pumped charge Qp 
in unit of 2e. An important measure of the accuracy of 
the pump is given by the deviation of the pumped charge 
from the quantized value 



In addition, since Cooper pair pumping is a coherent pro- 
cess, it is also interesting to test the accuracy of the opti- 
mization protocol by measuring the departure of the final 
quantum state from the desired one. To this end we also 
study the fidelity at the end of time evolution 

^^|((l|®(G(0)|)|*(T)}p 

which measures the overlap between the final state of the 
array and the ground state |G(0)}. In all plots we show 
the infidelity defined as X = 1 — JF. 



IV. OPTIMAL COOPER PAIR PUMPING 

In this section we present our results for the non- 
adiabatic optimal pump with if = 0, i. e. when there is 
no supercurrent. The goal is to achieve a pumped charge 
as close as possible to 2e (we remind that <p = is the 
case where the error is maximal without optimization). 



A. Optimal pumps 

Motivated by the adiabatic CcLSG, clS initial guess for 
control parameters we choose a circular path in (91,(72) 




FIG. 3: (Color on line) Top panel: Error in pumped charge £ 
and Bottom panel: Infidelity I after one cycle as a function 
of radius of circular path 50 in {qi, q2) space, for Ec /Ej — 10 
and 100. The vertical axis in both panels is in logarithmic 
scale. Results are obtained after using the quantum optimal 
control theory to modify the circular path (at least 200 iter- 
ations). T — SOnh/Ej and ip — 0. Numerical error, due to 
discretization of the time evolution, is at most of the order of 
10"^ 



space, around the degeneracy point 91 = 92 = 5 with 
radius qo described by sinusoidal gate voltages qi = 
i -f go cos(27rt/r -I- 6*0) and 92 = 5 + go sin(27rt/T -I- 6*0), 
where — 57r/4. The number of excess charges consid- 
ered for each box is n = 0, ±1, so that nine lowest charge 
states of the whole system are allowed to contribute to 
pumping. Since the initial paths with radii greater than 
1/2 might involve more states, we take go < 1/2. For 
each given radius, optimal control algorithm modifies 
the initial path such that the fidelity after one cy- 
cle approaches the value one. In the following we shall 
mostly focus on an experimentally relevant case and take 
Ec/Ej = 10. In this case in the adiabatic regime the 
accuracy of the pumped charge is always smaller than 
fifty percent^. As shown in Fig. [3l accurate pumping 
could be achieved with T = SOirh/Ej, which is already 
in the non-adiabatic regime, after at least 200 iterations. 
In the top panel of Fig. [3] we plot the error in pumped 
charge f as a function of radius go, while the bottom 
panel shows the infidelity X. Both error and infidelity 
are always less than 10~^. In Fig. [3] we show also some 
results for Ec/Ej = 100 (triangles), with error and infi- 
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FIG. 4: (Color on line) An example of gate voltages qi (top 
panel) and q2 (bottom panel) as a function of time, before 
and after optimization (200 iterations). The solid red lines 
correspond to initial circular path with radius qo — 0.275 
around the degeneracy point qi = q2 ~ 1/3 in {qi,q2) space. 
Dotted blue lines are optimized voltages leading to accurate 
pumping. T = SOnh/Ej, Ec/Ej = 10 and (p = 0. 
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FIG. 5: (Color on line) An example of Fourier transform of 
optimized gate voltage qi (blue curve in the top panel of figllll . 
Top panel shows most important harmonics as a function of 
frequency and few first harmonics are shown in bottom panel. 
qa = 0.275, r = ZQ-Kh/Ej, Ec/Ej = 10 and = 0. 



delity always smaller than 10^^. The difference between 
the accuracy of these two cases is due to the stronger 
violation of adiabatic condition, T 3> hEc/Ej, in the 
case Ec/Ej = 100 compared to Ec/Ej = 10. Notice, 
moreover, that the optimization procedure allows to con- 
trol directly the fidelity but only indirectly the value of 
pumped charge so that the behavior of T and f as a 
function of qq are not expected to be equal. 



B. Imperfections in the pulse shapes 

It is important to understand to which extent the 
pulses which optimize the pumping can be realized in 
experiments. Fig. [4] shows an example {qq — 0.275) of 
initial (red line) and optimized (blue dotted line) gate 
voltages qi (top panel) and q2 (bottom panel) after 200 
iterations. In the top panel of Fig. [5]we show the Fourier 



transform qi{uj) of the optimized gate voltage qi(t) plot- 
ted in Fig. m while in the bottom panel the same curve is 
plotted on a smaller range of frequencies. Note that the 
largest contribution occurs for lo = ±l/15Ej/h^ which is 
the frequency of the initial sinusoidal gate voltages. In 
a realistic situation, however, it is very difficult to imag- 
ine that all the details of the pulse, encoded in the high 
frequency components, can be reproduced faithfully. We 
accounted for imperfections in the pulse shape by intro- 
ducing a bandwidth parametrized by a high frequency 
cutoff uJcutoff- In Fig. [6] we show both the error and the 
infidelity as a function of the bandwidth for three differ- 
ent radii of initial circular path. No changes occur by de- 
creasing UJcutoff from 100 Ej/h until uJcutoff ~ 15Ej/h 
is reached and at this point a dramatic increase in both 
error and infidelity occurs which demonstrates the impor- 
tance of harmonics with frequencies close to this value. 
Another change in pumped charge and fidelity happens 
when UJcutoff ~ ^Ej/h is reached. Therefore the most 
important harmonics for optimizing the gate voltages lie 
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FIG. 6: (Color on line) The error £ in pumping one Cooper 
pair in each cycle (top panel) and the infidelity I of ending up 
to the ground state of the array at the end of cycle (bottom 
panel), in optimized non-adiabatic a function of the 

cutoff frequency uicutoff of gate voltages. Different curves are 
related to different radii of circular path in (91,52) space as 
the initial guess for control parameters in quantum optimal 
control theory. Ec/Ej = 10, T = SOtvU/Ej and ip ^ 0. 
Numerical error is at most of the order of 10"''. 



below uj = 20Ej/h and there is no need for higher fre- 
quencies in order to pump accurately. 



C. Effect of noise 

We finally discuss the effect of external noise. Since 
Cooper pair pumping is a coherent process, the presence 
of an external environment may be disruptive. In Joseph- 
son nanocircuits in the charge regime the dominant mech- 
anism of decoherence is 1// noise (see e.g. Ref,'^^). It is 
important to know whether the pumped charge is stable 
while optimized gate voltages are affected by noise. Al- 
though its understanding is far from complete, 1// noise 
is believed to originate from two-level fluctuators present 
in the substrate and/or in the insulation barrier. Sev- 
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FIG. 7: (Color on line) Error in pumped charge £ (top panel) 
and infidelity I (bottom panel) after one cycle as a function 
of radius of circular path qo in (91,52) space, for Ec/Ej = 
10 and 100, after applying the noise with power spectrum 
Sqf.{uj) = A/lo and strength A — 10~^ on optimized gate 
voltages. T — 30nh/Ej and (fi — 0. The vertical axis in both 
panels is in logarithmic scale. Numerical error is at most of 
the order of 10"^ 



eral theoretical works have recently studied the effect of 
1// noisc"^^. Following current approaches, we model the 
environment as a superposition of bistable classical fluc- 
tuators resulting in an additional random contribution 
6qk{t) to the gate charges after optimization. A distri- 
bution of switching rates 7 behaving as ^(7) oc I/7 in 
a range [7mm : Jmax] results in a noise power spectrum 
Sq^{uj) — {Sqk{t)Sqk{0))uj ~oj^^- We chose the switching 
rates such that the 1// part of the spectrum is centered 
around the typical frequency of the pump. Typically the 
1/ f region extends over two order of magnitudes in cu. 
We assume that fifty independent fluctuators are cou- 
pled weakly to the system and that the charge noise on 
the two separate gates are uncorrelated. Moreover we 
averaged the results over fifty different configurations of 
noise. In Fig.[7|we plot the error £ (top panel) and infi- 
delity I (bottom panel), as a function of qq, for a strength 
of noise A — 10^^ (^(w) — A/uj), which is a typical ex- 
perimental value. Remarkably, the pumped charge is vir- 
tually unaffected by noise: for Ec/Ej = 10 the error in 
the pumped charge and infidelity are still less than 10~^, 
which clearly shows the stability against noise. The case 
Ec/Ej = 100 is also stable. To complete the analysis. 
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FIG. 8: (Color on line) The error £ in pumping one Cooper 
pair in each cycle (top panel) and the infidelity X of ending up 
to the ground state of the array at the end of cycle (bottom 
panel), in optimized non-adiabatic function of the 

strength of the 1// noise A affecting gate voltages. Different 
curves are related to different radii of circular path in (gi, 52) 
space as the initial guess for control parameters in quantum 
optimal control theory. Ec/Ej = 10, T = SOivh/Ej and 
ip = 0. Numerical error is at most of the order of 10~^. 



we plot in Fig. [5] the error in pumped charge (top panel) 
as well as the infidelity (bottom panel) as a function of 
A for a few different values of go (note that both vertical 
and horizontal axes are in logarithmic scale). For all val- 
ues of radius go the accuracy of pumping and fidelity are 
larger than 90% even under the effect of noise with signif- 
icant strength A = 10~^, which is a good improvement 
considering that the accuracy is less than fifty percent in 
adiabatic regime without noise. More importantly, both 
£ and 2 are not increasing up to about A = lO^'^. 



D. Imperfections in the pump 

Up to now we have assumed the ideal situation in 
which the parameters of the pump are known, fn this sec- 
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FIG. 9: (Color on line) The error £ (top panel) and the infi- 
delity X (bottom panel), as a function of the random fraction 
xo added to the nominal Ec and Ej. Optimized pulses ob- 
tained for uniform array are used to calculate pumped charge 
and fidelity Ec/Ej = 10, T = SOttH/ Ej and ip ^ 0. Numer- 
ical error is at most of the order of 1Q~^. 



tion we address the effect of imperfections arising when 
such parameters are known only up to a given uncer- 
tainty from a uniform array. More precisely, we estimate 
how error in pumped charge and infidelity increase, as a 
function of the extent of the uncertainty, by evolving the 
system using the optimized pulse shapes calculated for 
the uniform array. To assess this issue we have imple- 
mented such an uncertainty on the parameters with two 
methods. In the first one we add to the nominal Ec and 
Ej, of a uniform array, a random fraction x taken in the 
range [— a;o,a;o]. In Fig. [9] we plot error and infidelity, re- 
spectively, averaged over 50 configurations as a function 
of xq for Ec ^ lOEj, for go = 0.300 and go = 0.475. In 
the second method we perturb the uniform Hamiltonian 
by adding to it a weighted random matrix. In Fig. [TU] 
we plot error and infidelity, respectively, as a function of 
the weight xq for Ec = WEj and for two values of go. 
Both methods show that no increase of infidelity and ep- 
silon is found up to a fraction (or a weight) of the order 
of 10~^. Finally, we emphasize that we have always as- 
sumed a uniform array simply for the sake of definiteness. 
Indeed, in the case of non-uniformity of the parameters, 
as long as they are known, the pulse optimization can 
always be performed. Of course, the shape of the final 
optimized pulses would be different, with respect to the 
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FIG. 10: (Color on line) The error £ (top panel) and the 
infidelity T (bottom panel), as a function of the weight xq of 
the random Hamiltonian added to the uniform Hamiltonian. 
Optimized pulses obtained from non-perturbed Hamiltonian 
are used to calculate pumped charge and fidelity. Ec/Ej = 
10, T = SOnh/Ej and ip = 0. Numerical error is at most of 
the order of 10~^. 



a conventional gate-based Cooper pair pump with op- 
timized pulse shapes. This pump operates in the non- 
adiabatic regime and, as we showed, can achieve a high 
accuracy. Numerical results after using the optimal con- 
trol theory were presented demonstrating that pumped 
charge is accurate and stable against the noise applied on 
gate voltages and uncertainty in charging and Josephson 
energy. Moreover all important harmonics contributing 
to optimized gate voltages have frequencies less than few 
Ej/h. The simple design of the pump, the short time 
scales of operation and the stability against noise are ad- 
vantages which make it worthy to think of the quantum 
optimal control theory as a powerful tool to achieve more 
realistic and still accurate pumping. As far as experi- 
ments are concerned, the limits on the bandwidth of the 
pulse generator are still demanding. 

In order to apply optimal control theory we enlarged 
the Hilbert space to take into account a passive detec- 
tor which acts as a counter. We assumed that initially 
the counter was set to zero implying that we supposed to 
disconnect the superconducting network from the elec- 
trodes. Moreover, we concentrated only on the case 
ip = 0. It would be important to find other methods 
for optimization which do not need to introduce counter 
so to explore also phase dependent errors. 



uniform situation, depending on the extent of the non- 
uniformity. 
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